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Abstract 

We present results on the QCD phase diagram for < vrT. Our simula- 
tions are performed with an imaginary chemical potential fi = ifij for which the 
fermion determinant is positive. On an 8 3 x 4 lattice with 2 flavors of staggered 
quarks, we map out the phase diagram and identify the (pseudo-)critical tem- 
perature T c (uj). For Hi/T < ir/3, this is an analytic function, whose Taylor 
expansion is found to converge rapidly, with truncation errors far smaller than 
statistical ones. The truncated series may then be continued to real fj,, yielding 
the corresponding phase diagram for hb ^ 500 MeV. This approach provides con- 
trol over systematics and avoids reweighting. We compare it with other recent 
work. 



1 Introduction 



A prime goal of heavy ion collision experiments at SPS, LHC (CERN) and RHIC 
(Brookhaven) is to probe the transition^ from hadronic matter to a quark gluon plasma 
at high temperatures and small baryon density Because QCD is strongly interacting, 
the only first principles method to predict the phase diagram is by means of lattice 
simulations. On the other hand, at finite densities the fermion determinant is complex, 
thus prohibiting standard Monte Carlo sampling of the path integral in what is known 
as the 'sign problem'. Overviews with references to various studies of this problem may 
be found in 

Recently, a first (/i, T) phase diagram for 2+1-flavor QCD has been presented by 
using a two-dimensional generalization || of the Glasgow reweighting method @]. For 
such methods cancellations in the reweighting factor occur generically at large volumes 
and/or chemical potentials, and it remains to be seen how well this approach allows 
for infinite volume and continuum extrapolations. Another work |J attempts to im- 
prove on this aspect by expanding the product of reweighting factor and observable 
in powers of (/i/T), and computing the coefficients through second order. The calcu- 
lation is then equivalent to computing susceptibilities at fi = ||, which is possible 
on any volume. The /i/T-range of applicability is limited by either the /i/T- value of 
the expected critical point marking the onset of a phase transition, where suscepti- 
bilities become non-analytic, or the radius of convergence of the expansion, whichever 
is lower. However, a priori nothing is known about the convergence of the expan- 
sion or the error introduced by its truncation. Moreover, a problem common to all 
reweighting approaches is that so far it has been impossible to assess the overlap of 
reweighted ensembles with the original. Hence, the results obtained in 0, |5| rest on as 
yet uncontrolled approximations. 

In this paper we attempt to develop a controllable alternative approach, avoiding 
reweighting in /i altogether. This is achieved by simulating with imaginary chemi- 
cal potential, where there is no sign problem and hence no need for reweighting. We 
employ the idea |7]] that in this case one may fit the non-perturbative data for an ob- 
servable by truncated Taylor series in /i/T, thus keeping full control over the associated 
systematic error. In the parameter range where this is possible, and in the absence of 
non-analyticities, the series may be analytically continued to real values of /i. Non- 
perturbative evidence for the viability of this approach has been given for screening 

^Jn this paper the term "transition" is used to label a change in dynamics, whether it proceeds by 
a smooth crossover or by a first or second order phase transition. For the last two, we use explicitly 
the term "phase transition". 



1 



masses in the deconfined phase ||, which can be simulated sucessfully for real and 
imaginary // [R] in the framework of dimensionally reduced QCD [pTD| . Here we extend 
this approach to the critical line T c {fi). This is possible by noting that, while giving the 
location of singularities in thermodynamic limit, the critical line itself is an analytic 
function. After continuing it from imaginary to real \x for various volumes, a finite 
volume scaling analysis of its location should also reveal the nature of the transition. 
Since there is no sign problem, there are no limitations on the volumes that can be 
simulated either. In particular, it should then be possible to identify the location of 
the critical endpoint of the phase transition which is expected on general theoretical 
grounds [11], and for which reweighted numerical results have been presented in [[|. 

Our method is limited by a maximal imaginary part of the chemical potential, /if = 
7rT/3, above which the system tunnels into an unphysical Z(3) sector. Along the 
transition line, this corresponds to an accessible range of real /i£<500 MeV, which 
well covers densities in heavy ion experiments [1^]. We present numerical results for 
the location of the critical line, and find that the truncation error of the Taylor series is 
far smaller than statistical errors in the whole accessible range. A finite volume scaling 
analysis is beyond the scope of our current computational resources and postponed 
to later work. Nevertheless, our results indicate that the part of the phase diagram 
relevant for heavy ion collisions can be reliably obtained from simulations of imaginary 
chemical potentials. 

In Sec. |2] we recall some general features of QCD at imaginary chemical potential, 
which will be useful in the sequel. The analyticity of the critical line is derived in Sec. [| 
Numerical results for two flavor QCD with Wilson gauge action and Kogut-Susskind 
fermions are presented in Sec. §, while Sec. |^ discusses our approach in comparison 
with previous ones. In Sec. |6| we draw our conclusions. 



2 QCD at real and imaginary chemical potential 

While this section contains no new results, we would like to discuss some qualitative 
features in detail as they play an important role for the following. We denote by \i the 
chemical potential for quark number Q, and by \xb the chemical potential for baryon 
number B = Q/3, i.e. fi = /is/3. The grand canonical partition function is compactly 
written as 

Z(V,fi,T) = Tr (e-(6-»QWry (1 ) 

An expansion for small chemical potentials proceeds in terms of the dimensionless 
parameter Jjl = fi/T, and hence our interest is in Z(V,p,,T). 
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The presence of a chemical potential term in the action breaks the separate invariance 
under Euclidean time reflection and charge conjugation. In particular, a time reflection 
can be compensated by \l — > —\l and vice versa, so that 

Z{p) = Z(-fl). (2) 

As a consequence, observables symmetric under time reflection are even functions of Ji. 
Let hi G R denote the real and imaginary parts of \x\ 

H = fi R + i^j, (3) 

and similarly for /2. The qualitative features of SU(N) QCD at imaginary chemical po- 
tential, lir = 0, were studied a long time ago With dynamical fermions the Z(N) 
symmetry of the pure gauge theory is explicitly broken. However, in the presence of a 
complex chemical potential, a Z(N) transformation of the fermion fields is equivalent 
to a shift in lii, thus leading to a new symmetry: the partition function is periodic in 
fjij with period 2irT/N @, 

Z(fi R ,Ji I ) = Z(jl R ,jl I + 2ir/N). (4) 

The different Z(N) sectors are characterized by the phase <p of the Polyakov loop, 

(P(x)) = \(P(x))\&. (5) 

At high temperatures one may use the perturbative effective potential V(ip) to obtain 
the qualitative behaviour of the theory. In pure gauge theory (P(x)) = for T < T c 
and (P(x)) for T > T c . In the deconfined phase there are N degenerate minima at 
ip = 2irk/N, k = 1, . . . N, separated by potential barriers, and the transition between 
them is of first order |14|| . For Nf > 0,fi = 0, explicit Z(iV)-breaking results in an 
expectation value for the Polyakov loop even in the confinement phase. The expectation 
value is real, tp = 0, corresponding to the true vacuum of the theory, while the other 
sectors are metastable [|15|]. This is true both for T > T c and T < T c . When a chemical 
potential is switched on |13[ [16|, 17[], its real part stabilizes this situation by raising the 



energy of the other Z(N) vacua, while the imaginary part has the opposite effect of 
lowering one of them relative to the (p = vacuum. Hence, once /// exceeds some critical 
value fi c j, a phase transition to a non-trivial Z(N) sector occurs. Non-perturbatively 
the same statement follows from the aforementioned equivalence of certain /^/-shifts 
and Z(N) rotations. The combined symmetries Eqs. (0),(ffl imply that this transition 
is periodically repeated for larger values of Moreover, its exact critical values are 
dictated to be [|13| 



N V 2 
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On the other hand, for purely real /x this transition never occurs. This situation is 
sketched qualitatively for SU(3) in Fig. [I]. The order of the Z{3) transitions depends 



71 



tt/3 



Figure 1: Location of Z(3) transitions, Eq. Analytic continuation is limited to 
the region within the arc. 

on temperature. A perturbative calculation predicts it to be of first order in the 
deconfined phase, whereas at low temperature a strong coupling analysis suggests it 



to be a smooth crossover |13|]. We shall present numerical results in support of this 
picture in Sec. f|. 

Note that the first of these Z(3) transitions at pfj = tt/3 limits the useful information 
to be obtained from imaginary chemical potentials: Eqs. (§),([§) imply 



Z(fii = tt/3 + A/ij) = Z(jii = 7r/3 — A//j) 



(7) 



so that all observables are symmetric about the transition at £l c j. In other words, all 
expectation values for pj > fij are simply copies of those for JJlj < p,j, determined by 
reflection symmetry and periodicity. Moreover, here we are interested in analytically 
continuing results obtained on the imaginary axis, and this is only possible for \p,\ < tt/3 
||, represented by the arc of circle in Fig. [[]. Furthermore, the observable to be 
continued has to be free of additional singularities in this region. In the next section 
we will show that this is the case for the critical coupling for deconfinement, and hence 
the associated temperature T c (/i). Along the deconfinement transition the latter is in 
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the range T ~ 160 — 170 MeV, which means that our approach is restricted to baryon 
chemical potentials //# < 500 MeV. 



3 Analyticity of the (pseudo-) critical line 

On a L 3 x N t lattice, the rescaled chemical potential takes the form p, = N t afi, which 
up to the constant factor N t is equivalent to the chemical potential in lattice units. 
On the lattice one then explores the phase diagram in the (/3, /i)-plane, which in the 
continuum limit has to be converted to a (T, /i)-diagram. 

The location and nature of a phase transition may be determined from the finite 
volume scaling of susceptibilities of dimensionless operators (V = L 3 ), 

X = VN t ((0-(0}) 2 ), = ^-J20(x). (8) 

In pratice we shall use as 0(x) the plaquette, the chiral condensate and the modulus 
of the Polyakov loop. In a transition region, fluctuations are strong and susceptibilities 
display a peak, Xmax = x{Pc,Pc), whose location determines the critical parameters. 
In a finite volume, the susceptibility is always an analytic function of the parameters 
of the theory, even in the presence of a phase transition. The latter reveals itself by a 
divergence of Xmax in the infinite volume limit, whereas Xmax stays finite in the case of 
a crossover. The order of the transition is determined by the rate of divergence: for a 
first order transition Xmax oc V, whereas for a second order phase transition, Xmax oc V p 
with a critical exponent p = ^/du < 1 |fTB| . Alternatively, one may consider the finite 
size scaling of the critical coupling (3 C (V) itself, which attains its infinite volume limit 
as 

(P c (V)-[3 c (oo))~V-°, (9) 

where a = 1 for a first order phase transition, o = Xjdv < 1 for a second order phase 
transition, and a = for a crossover. (For notation and numerical application to SU(3) 



pure gauge theory, see e. g. jn| and references therein). 

At zero density, the deconfinement transition is signalled by a peak Xmax = x(0,/5 c ), 
defining the critical coupling through the equations 



dx 
dp 



d 2 x 





dp 2 



< 0. (10) 



When the chemical potential is switched on, this peak extends as a "mountain ridge" 
into the (/?, p) -plane. Starting at /3 C (0), the ridge is defined by maximizing x while 
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moving into the (/?, /i)-plane, i.e. by moving along the direction of the gradient Vx(j£t, /3) 
(or the direction with the largest second derivative in the case of vanishing gradient). 
It has a vanishing first and negative second derivative in the direction perpendicular 
to the gradient. 

For the numerical analysis, however, it is more convenient to consider maxima in the 
directions of (3 and /2, respectively. The two pairs of conditions 

{fyx = 0,^x<0}, {0 A x = O,9jx<o} (ii) 

then define two separate lines bounding the ridge from below and above, as indicated 
in Fig. |2] (left), where the ordering of the lines depends on the sign of the gradient. 
For large enough volumes one expects a positive gradient while increasing \i towards a 
critical point, and the first condition denotes the lower line. When the volume is taken 



T 



finite V 




T 



infinite V 



Figure 2: Schematic location of the lines defined by Eq. (fH[). In the infinite volume 
limit (right), the lines merge for a phase transition, but stay separate for a crossover, 
the bifurcation marking the end point of the phase transition. 



to infinity, the mountain ridge changes its shape: where there is a phase transition, 
it becomes singular and infinitely narrow, i.e. the two lines merge into one. Where 
there is a crossover, the ridge saturates at finite height with finite curvature, and 
the two lines remain separate (as in Fig. |^ (right) at small /i). In this case there is 
no phase transition and no uniquely defined pseudo-critical line. Nevertheless, the 
physical properties change rapidly in the region of the line defined by dp\ = 0, as is 
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well known from the zero density case. Moreover, for small p this line is the one closer 
to the ridge. Thus, for every p we may locate the thermal transition by the conditions 
specified in Eq. (|T0|), which then represent our implicit definition of the critical line 

Pc(fi). 

Let us now consider a finite volume. Expanding the susceptibility around ft c (0) and 
using xifii ft) — x(— fi->ft), h takes the form 

X (p,ft) = E c nm ((3-(3 c (0)) n fi 2m - (12) 

n,m=0 

Because x(/Z, (5) is analytic, its derivatives are as well. In the neighbourhood of the 
critical line we furthermore have d^x 7^ 0- The implicit function theoremf] then tells 
us that the implicitly defined ft c {fi) is analytic in p. Furthermore, its first derivative is 
obtained from the chain rule for partial differentiation, 



dp dpd(3 [dp 2 J 



(13) 



Since x is even in fi, the same follows from the last equation for /3 C , which therefore 
has a Taylor expansion 

ftc(p) = E a nf n = E Cn(a/i) 2n , (14) 
n=0 n=0 

where we have absorbed the temporal lattice length into the coefficients in the second 
equation, c n = a n N 2n . 

All of these considerations are unchanged if we consider a purely imaginary potential, 
p, = ipLi. In this case we may define a corresponding real function xifiii ft) with a Taylor 
expansion in ft,p^ and coefficients c nm . Clearly, this is just the analytic continuation 
x(P>i,ft) — x(ip>iift)i with c nm = c nm (— l) m . Its maximum in the /3-direction is again 
defining a function /3 c (/2j), which by the implicit function theorem is analytic. The 
analogue of Eq. (]TJ) then reads 

m d 2 x (d'xY 1 . d\ (d'xY 1 _.dftc . , 

d\ki dfnd(3\dp 2 ) l dfxdp\dl3 2 ) 1 dp, ' 1 ' 

Since this equation is true for every pi < pj we have thus established ft c (pi) — ftc(P = 
ipi), i.e. the critical line of the deconfinement transition at imaginary chemical potential 
is simply the analytic continuation of the critical line at real chemical potential. 

2 See e.g. |2(J. We thank F. Wilczek for reminding us of this theorem in this context. 
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Next, we need to discuss the infinite volume limit. As remarked above, for every /2 
the critical coupling will approach its thermodynamic limit in a fashion characteristic 
of the order of the transition. However, /3 c (oo, Ji) is only shifted from (3 C (V, p), without 
developing any singularities in the thermodynamic limit. Hence, the critical line defined 
by Eq. ( pf) marks the singularities in the partition function, which it smoothly connects 
to a pseudo-critical line in the crossover region. It thus remains itself an analytic 
function of ft in the thermodynamic limit for all ft. 

Our strategy is now to first compute the critical line for imaginary chemical potential 
and check the convergence of its Taylor expansion. As we shall see, convergence is 
surprisingly fast for the whole range of chemical potentials accessible to this method. 
Analytic continuation then reduces to simply flipping the sign of the appropriate terms. 
Finally, the infinite volume limit has to be taken from the continued results at real 
/i. The order of the transition can then be determined in a (V, /2)-range where the 
truncation error of the series is smaller than finite size scaling effects. 



4 Numerical results for two light flavors 

In order to test the feasibility of our approach, we consider QCD with two flavors of 
staggered fermions with bare mass am = 0.025 and = on a 8 3 x 4 lattice. We use 



the R-algorithm |2J] with a step size St = 0.02, sufficiently small for the systematic 
errors 0(St 2 ) to be negligible compared to our statistical errors. For 8 values of (3 
spanning the relevant temperature regime and 6 values of afij e {0,7r/12}, we have 
accumulated 2000 — 6000 unit-length trajectories each, measuring the gauge action, the 
Polyakov loop and the chiral condensate after each trajectory. For one parameter set 
we have performed additional simulations on 6 3 x 4. To determine the pseudo-critical 
value p c (afii), we use the Ferrenberg-Swendsen reweighting method 



4.1 The Z(3) transition for imaginary fi 

First, it is instructive to investigate the Z(3) transition discussed in Sec. [| The pre- 
dicted critical value for the chemical potential is at fifj = tt/3, which on our N t = 4 
lattice corresponds to a/ij = 7r/12 = 0.262. For this value Fig. |3] shows histograms of 
the phase <p of the Wilson loop, obtained for various lattice couplings (temperatures). 
We find the system tunneling between the sectors with ip = and if = — 27r/3. 

In accord with the predictions made in we observe a smooth distribution consis- 
tent with a crossover at low temperatures, and a pronounced two-state signal indicating 
a first order phase transition at high temperatures. The critical coupling for the onset 
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Figure 3: Probability distribution of the phase of the Polyakov loop for the critical value 
ajjfj = n/12. 
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Figure 4: The (3- and fij- dependence of the average phase of the Polyakov loop. 

of the Z(3) transition is visibly between 5.313 < (3 C < 5.34. 

For low temperatures, where the Z(3) transition is continuous, there is frequent 
tunneling between the Z{3) sectors. It then follows that also for a\ii < (afii) c the 
ensemble average consists of a mixture of the two sectors, where the ip = — 27r/3 
admixture grows with a/i/. This is clearly visible in Fig. |], where the average phase 
of the Polyakov loop gradually moves from (p(afii = 0) = to ip((afii) c ) = —n/3, the 
latter value corresponding to the average between the two phases ip = 0, — 2n/3. For 
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0.5 1 

3Hj/(27tT) 

Figure 5: Schematic phase diagram for imaginary chemical potential: The vertical line 
marks the Z{3) transition (this section), the curved lines the deconfinement transition 
(next section). The solid line indicates a first order transition, while the nature of the 
dotted lines is not yet determined. The diagram is periodically repeated for larger values 
of fr. 

a/if > (a/i/) c , the phase gets dominated by the vacuum with tp = —2it/3. 

At high temperatures, on the other hand, the Z{3) sectors are separated by a di- 
verging potential barrier, and tunneling is more suppressed the higher the temperature. 
As a result the Monte Carlo ensemble stays longer in one sector, and its equilibrium 
is dominated by the lowest vacuum. Thus, as shown in Fig. [|, for a/zj < (a/i/) c the 
system always settles in the trivial vacuum with zero phase, whereas for a\i > (a/ii) c 
it settles in the vacuum sector with tp = —2n/3. The dividing line between those 
situations is ip{{afii) c ) = —tt/3, which is independent of temperature by symmetry. In 
Fig. B, one observes that indeed p((afj,i) c ) = —tt/3 for all /3's, although errors become 
large at large (3. This reflects the difficulty of tunneling at high temperatures, and 
implies longer and longer Monte Carlo time for the system to equilibrate. Hence, at 
high temperatures great attention must be paid to Z{3) ergodicity and its possible 
violation. 

From the susceptibilities x(P) a t (a/i/) c = 7r/12, only one peak is apparent. This 
means that, within our accuracy, at (a/i/) c the critical temperature at which the Z(3) 
first-order transition line starts is the same as the critical temperature for the deconfine- 
ment transition determined below. Schematically the (T, /j/) phase diagram then looks 
as in Fig. [5], with a special point where three critical lines meet. This point might 



10 



a (3=5.30 
□ (3=5-32 
o (3=5.35 



Oh 

x> 



Figure 6: The chiral condensate ipip as a function of a\ii. 

be tricritical if deconfinement corresponds to a genuine phase transition, or critical 
otherwise. We take a first look at this issue below. 

Once ergodicity is ensured, the effects of the Z(3) transition are also visible in other 
observables like, e.g., the chiral condensate, as shown in Fig. || Note the symmetry 
of the observable around afij according to Eq. (0). For low (3 (viz. temperatures) the 
observable is continuous at the critical chemical potential, while a cusp indicating a 
discontinuous transition develops for large f3. 

Let us now try to determine the endpoint (f3 c , a^i = vr/12) of the first-order transition 




Figure 7: Ratio of of cumulants (<p ) / {<f ) , for a\ii = tt/12 
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line with better accuracy. We can do so by monitoring the plaquette susceptibility as 
a function of (3. This procedure gives us f3 c = 5.325(5). For improved accuracy, and to 
check our control of the delicate ergodicity problems mentioned above, we performed 
additional simulations on a 6 3 x 4 lattice at a/// = tt/12. Since ip = — vr/3, one can 
form a ratio of cumulants (ip 4 ) / (ip> 2 ) 2 , where <p = (p + ~, and estimate the critical 
point (3 C from the crossing of the cumulant ratios for the two available volumes. This 
procedure, illustrated in Fig. [7], yields a consistent value f3 c = 5.321(7). Note also that 
the crossing occurs at a value of the cumulant ratio 1.8(2) consistent with that of the 
3d Ising universality class 1.604(1) |23[. This indicates as the simplest interpretation of 
the phase diagram: the first-order Z(3) transition line ends with a second-order critical 
point in the Ising class. This picture will likely depend on the number of quark flavors 
and on their masses. 
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4.2 The deconfinement transition for imaginary /i 

In order to map out the deconfinement transition line in the (/?, /2/)-plane, we have 
measured the susceptibility of the plaquette, of the modulus of the Polyakov loop, and 
the disconnected part of a stochastic estimator for the chiral condensate susceptibility. 
All three observables yield consistent results for /3 c (a//j), as shown in Fig. ||. We have 
chosen to display only three a/ij-values for clarity of the plots. Note that the critical 
lattice coupling is growing with increasing a///, and thus the critical temperature is 
an increasing function of /J? as well. This is in contrast to the situation at real 
where T c (/x) is a decreasing function. This qualitative behaviour is in accord with the 
one found for screening lengths in units of inverse temperature, which are increasing 




P 

Figure 8: Susceptibilities for the plaquette, the chiral condensate and the absolute value 
of the Polyakov loop. 
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Figure 9: Location of the pseudo-critical deconfinement line as determined from pla- 
quette susceptibilities. The lines show the fits given in Table [J. 
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Table 1: The first coefficients of the Taylor expansion of the critical coupling, Eq. flF${). 
Data and fits for /3 C from Xpiaq are shown in Fig. |^. 

(decreasing) functions of \Lj (hr) |§. (An increasing screening length means more 
energy is required to separate charges before they are screened, and hence an increasing 
deconfinement temperature). 

Since the plaquette susceptibility is the one that is most accurately determined, it 
is our choice for determining the critical line. The values of j3 c corresponding to six 
values of afij are shown in Fig. [5| Our next task is to determine the reliability of 
the Taylor expansion for the critical line, Eq. (|I4]) . Results of fitting the data to a 
polynomial of degree one and two in (a/i/) 2 , respectively, are shown in Table [I| and 
Fig. p. Excellent fits are obtained in both cases with compatible coefficients c^i- 
The quartic correction appears negligible even near a/ij = (a/i/) c , leaving c 2 to be 
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consistent with zero within this range. This is again in accord with the analogous 
finding for screening masses, which are equally well described by the leading (a///) 2 
term in that range || . Statistically consistent results are obtained when f3 c is extracted 
from Polyakov loop instead of plaquette susceptibilites, as also shown in Table [l]. We 
then conclude that the (pseudo-) critical line for imaginary chemical potential ft = ifti, 
fti < 7r/3, is well described by 

/3 c (a/i/) = c + ci(a/i 7 ) 2 , (16) 
with Co, ci as per Table [l] (top line). 



4.3 The deconfinement transition in physical units for real /i 

With the results of the previous section, it is now trivial to obtain (3 c (ft) for real ft, by 
continuing ftj — > —ifti. As a comment on systematics, let us add that the difference 
between the linear and quadratic fits in (a/z/) 2 from Table [I] gets slightly amplified by 
analytic continuation (only the sign of one term is flipped). This reflects the fact that 
convergence properties in the real and imaginary directions may be different in general. 
Such an effect will be large for an observable with pronounced structures, and small on 
a smooth observable like ours. In any case, the effect can be monitored, and the error 
band of the first fit is even after continuation completely contained within the error 
band of the second fit, which can hence be dropped as statistically insignificant. 

It is instructive to translate our result into physical units, in order to illustrate the 
phenomenological relevance of our approach. However, we caution that at this stage 
such a conversion is merely illustrative, with quantitative numbers still afflicted by 
various systematic errors: while we have monitored that the truncation error of the 
Taylor series is negligible, there might still be sizeable corrections from finite volume 
and lattice spacing effects, and we have not investigated the quark mass dependence 
of our results. In view of this we content ourselves with the perturbative two-loop 
expression for the lattice QCD scale in the presence of two massless flavors, 

-bi/2b 2 

a((3)A L =1^1 exp(-/3/12& ), 



(17) 




2„\ , / 1 VA.„„ 38 

T c (/i) a{/3 c {fj))A L 



T c (0) a(/? c (0)A L 

while the scale is set by the critical temperature T c (/i = 0) = 173(8)MeV in the chiral 
limit for staggered fermions |24[ (which is expected to still contain finite lattice spacing 
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Figure 10: Location of the deconfinement transition corresponding to the first fit in 
Table [J. The error bar gives the uncertainty in T c (0) used to set the scale, the dotted 
lines reflect the error on ci from Table [/]. 



errors). This leads to the result for the (pseudo-) critical line as shown in Fig. [10]. As 
remarked before, a detailed finite volume study is beyond the scope of the present work 
but would reveal the order of the transition along the line. 



5 Comparison with other methods 

Before concluding, we would like to make some comments comparing our method with 
previous studies of the finite \x deconfinement transition [0, |5|. The main disadvantage 
of our approach is its limitation to the range \fi\/T < ir/3, as discussed in Sec. 0. On 
the other hand, within this range we have full control over the necessary approximation, 
i.e. the truncation of the Taylor series, which we find to pose no limitation at all within 
the accessible range. 

We would like to illustrate the importance of this aspect by a comparison with the 
reweighting method 0, 0. In Ref. 0, the overlap of the reweighted with the full 
ensemble is tested for imaginary \i. The configurations used for reweighting from the 
/i = ensemble are all in the real Z(3) sector, so that reweighting is by construction 
insensitive to the Z(3) transition. This is clearly seen in Fig. 1 of Ref. ||, where the 
reweighted chiral condensate is found to be a monotonically increasing function of (a/xj) 
in the range < a\ij < 0.32. Such a result is inconsistent with the symmetries of the 
imaginary-/x theory. Ref. || works at the same lattice spacing as we do, and therefore 
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the critical chemical potential for Z(3) tunneling is at (a/j,j) c = it/ 12. According to 
Eq. (0), the chiral condensate has to be symmetric about this point and decrease for 
larger a/i/, as it does in our data, Fig. [| Not surprisingly, reweighting completely 
misses this qualitative behavior|. In the case of real //, there is no Z(3) transition. 
However, the relevance of complex Z{3) sectors was studied in [[[7]]. It was argued 
there that constraining the Polyakov loop to be real is a reasonable approximation 
below and above the deconfmement transition, but a poor one near it. In any case, the 
discussion illustrates a potential danger of any reweighting method: it may completely 
miss a qualitative change in the physics without announcing its failure, which only 
shows up by comparison with results obtained without reweighting. 

Ref. studied Nf = 2 + 1 flavors, with light quark mass values and lattice spacing 
as ours, whereas the bare strange quark mass was chosen eight times heavier than the 
light ones. Keeping this difference in mind, we find the location of our continuum 
transition line consistent with theirs over the whole range we have studied. 

As mentioned in the introduction, another recent work || employs a Taylor expansion 
of the reweighted path integral. However, in this case only susceptibilities are measured, 
so that one has no control over the higher order terms. Let us recall that the Taylor 
expansion is in fi/T, which becomes larger than one already at fi ~ 170 MeV. A priori 
there is no way of knowing how large the convergence radius of such an expansion is. 
Simulating at imaginary \x provides a framework for obtaining quantitative information 
about the convergence. As it turns out, our results non-perturbatively endorse the 
use of the leading Taylor coefficient in T c (/x) up to \i ~ 170MeV. Ref. || works at 
quark masses four to ten times larger than ours, but finds only very weak quark mass 
dependence. Keeping this difference in mind, we find our result for the slope c\ of 
/? c ((a/i) 2 ) to be consistent with theirs. 



6 Conclusions 

The past year has seen progress in the numerical study of finite density QCD. How- 
ever, this progress is essentially based on adopting a more pragmatic viewpoint. The 
net quark density in heavy-ion collisions is small. Therefore, instead of solving the 
fundamental difficulty posed by the sign problem, it is very useful already to explore 
the regime of small chemical potential, where the sign problem is moderate and various 



3 Fig. 1 of Ref. H also shows the chiral condensate rising monotonically in the full ensemble. A 
possible reason for this strange behavior might be insufficient thermalization, cf. the discussion in 



Sec. 4.1. Indeed, at yet larger ///, the authors of Ref. H] find, as expected, clear disagreement between 



reweighted and full ensemble results 



17 



approaches appear possible. 

We have studied QCD with two light flavors of staggered fermions in the presence 
of an imaginary chemical potential /i = i/i/. In this case the partition function is 
periodic in fij and, in addition to the deconfinement transition, Z(3) transitions occur 
at (hj/T) c = 27r(n+l/2)/3. We mapped out the phase diagram in the (T, /z^-plane and 
found that the Z(3) transitions are of first order in the deconfined phase and continuous 
in the confined phase. For (fij/T) < n/3, the critical temperature for deconfinement is 
an analytic function, irrespective of the order of the transition. Within this range its 
Taylor expansion is found to rapidly converge and to be well described by the leading 
term. 

Other recent work uses reweighting of simulations performed at fi = 0, and at one 
or more temperatures. In contrast, we vary two parameters in our simulations: the 
temperature as usual, and also the imaginary chemical potential fij. Having at our 
disposal two-dimensional information allows us a reasonable control over systematic 
errors. 

We have used this property to determine the transition line T c (/z) for real fi through 
analytic continuation. While we cannot, on a single volume, find out whether the 
transition is a crossover or a genuine phase transition at each value of /z, a finite-size 
scaling study should elucidate this issue, provided that finite size effects are larger 
than the truncation error. Our central result is shown in Fig. [It]. It is consistent with 
Refs. 0, |5[. The transition line is well represented by the equation 



while the next-order term 0((hb/T) a ) is statistically insignificant up to fiB ~ 500 MeV. 
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